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ABSTRACT 

We present the previously unanalysed high quality XMM-Newton spectrum of an ultralumi- 
nous X-ray source candidate in NGC 45 17. As with other high quality ULX spectra, a down- 
turn in the spectrum is observed at ^6 keV. Both of the recent disc reflection and Compton- 
isation interpretations of this feature are applied, in order to present a direct comparison, and 
are found to provide statistically equivalent representations of the current data. We find that 
the reflection model requires the accretion disc to have a highly super-solar iron abundance, 
while the Comptonisation model requires low temperature Comptonising electrons, and for 
the corona to be optically thick. These physical requirements are discussed in detail, and 
physically motivated scenarios are highlighted in which each model can be considered a vi- 
able explanation for the observed emission. By extending our consideration of these two inter- 
pretations to high energies, we demonstrate that observations of ULXs at energies > 10 kcV 
should be extremely useful when attempting to distinguish between them. With current in- 
strumentation, it is only viable to perform these observations for M 82 X-1, but future high 
angular resolution hard X-ray imaging spectrometers, such as the Hard X-ray Imaging System 
due to fly on Astro-H, should go a long way to resolving this issue. 
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1 INTRODUCTION 

Ultraluminous X-ray Sources (ULXs), discove red with the Ein- 
stein observatory in the 1980s l lFabbiandlT98^) . are extra-nuclear 
point sources observed to be more luminous in X-rays than the 
Eddington luminosity for a stellar mass (~10Mq) black hole, 
i.e. Lx > 10'^^ erg s^^, a combination which has led to extended 
debate over the nature of these sources. The extremity of the ob- 
served luminosities implies that these sources are most likely to be 
accreting black holes of some kind. 

There are currently three main theories to explain the apparent 
X-ray luminosities of ULXs. The first is that these sources are sys- 
tems similar in nature to the less luminous X-ray binaries (XRBs) 
frequently observed, but in which the mass of the black hole is 
greater than the stellar- mass black holes that XRBs are expected to 
harbour. However, they must also be less massive than the super- 
massive black holes associated with active galactic nuclei (Mbh> 
10^"'' Mq; SMBH), as dynamical friction dictates that such ob- 
jects would sink to th e centres of their resp ective galaxies well 
within a Hubble Time jMiller & Colbertll2004) . In this case, ULXs 
are interpreted as interm ediate mass black holes with 10^ M© < 
Mbh < 10^ Mq (IMBH: IColbert & Mushotzkvll999h . The second 
is that ULXs are actually the same systems as regular XRBs, but ob- 
served in a different accretion state in which the black hole is able to 
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radiate at super-Eddington rates (for a rec ent review of the standard 
accre tion states observed in XRBs see iRemillard & McClintockl 
l2006l) . and a number of me thods have been proposed by which this 
might be possible (see e.g. IPoutanen et al ] |2007l , lFinke & Bottched 

The third again involves stellar-ma ss XRBs, but in t his case 
the observed emission is not isotropic jKing et aklllOOll) . ULXs 
may be such sources with which we have a favourable orienta- 
tion, and as such artificially high luminosities are calculated for 
these sources based on the assumption of isotropic emission. Vio- 
lation of the Eddington limit may be avoided by assuming a suit- 
ably high level of anisotropy. Such anisotropy could be due to 
geometrically thick accretion discs subtending a larger solid an- 
gle, funnelling the X-rays produced in the inner regions. Such 
'slim' discs ma y also allow a compact o bject to radiate at super- 
Eddington rates jAbramowicz et aljl980l) . so super-Eddington and 
anisotropic emission may be intrinsically linked. Alternatively 
anisotropic emission could manifest itself in the form of pencil- 
beam relativistic jets si milar to those seen in other accreting sources 
( iRevnolds et al.lll997h . However, a number of ULXs have been 
observed to be embedded within (roughly) spherically symmet- 
ric phot oionised emission nebulae inconsistent with narrow beam- 
ing (e .g. IPakull & Mirionill2003l iKaaret et al.ll2004 iBerghea et d] 
I2OI0I) . This suggests highly anisotropic emission by itself seems 
unlikely to be able to explain ULXs as a class, although it cannot 
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Figure 1. XMM-Newton X-ray (top) and DSS optical (bottom) images of 
NGC 45 17. The elliptical region shows the RC3 D25 isophote, and the cross 
and diamond indicate the positions of the nucleus and NGC 45 17 ULXl 
respectively as the position eiTor circles for the X-ray detections are too 
small to show. 



be mled out in individual sourc es. For a more detailed review on 
the possible nature of ULXs see iRobertsI ( |2007|) . 

Observationally it has been difficult to distinguish between 
the IMBH and stellar mass interpretations. ULXs are extragalac- 
tic sources hence detailed information on optical counterparts is 
still relatively rare, although this is being addressed, and orbital 
periods/parameters of the assumed binary system are even rarer 
still; attempts at dynamical measurements of the black hole mass 
are only now becoming possible for a small number of sources 
with the use of world leading observational facilities. Current 
examples of these a ttempts are limited to perio dicity claims in 
M 82 X-1 jKaaret & F eng 20 j), NGC 13 13 X-2 l lLiu et alj|2009l) 
and NGC 5408 X-1 dStrohmaveJ |2009|) . although in the latter 
case [Foster et al. suggest the periodicity might be super- 
orbital, related to jet precession. Much recent work has been fo- 
cused on successfully modelling the X-ray spectral and timing 
properties of ULXs, as both the inner temperature of the accre- 
tion disc and characteristic timescales are mass dependent, pro- 
portional to Mp,„^"^ and Mjj^ respectively l lMakishima et al.l200(l 
IVaughan et alj2003h . Early results pointed towards an IMBH inter- 
pretation with modelled disc temperatures being cool er and char- 
acteri s tic timescales being long er than XRBs (see e.g. iMiller et al.l 
1 20041 : IStrohmaver et alj |2007|) . However, a systematic study of 
the highest quality ULX X-ray spectra by IStobbart et"aLl | |2006|) 
showed the majority displayed a break/turnover at high energies 
(> 3keV). Such breaks are not commonly seen in the spectra of 
XRBs^ 

[Gladstone et al. I j2009l) argue that this turnover could be due 
to Comptonisation from an optically thick corona, which would 
shroud the inner regions of the accretion disc and artificially lower 
the inner temperatures obtained from simple spectral modelling. 
Under this assumption, the 'true' inner disc temperatures recov- 
ered are akin to those observed in XRBs. Most ULXs would then 
represent a new, high Eddington fraction accretion state for these 
ste llar-mass black holes. An alterna tive explanation is proposed 
bv lCaballero-Garcia & FabiaiJ bOlOl) who note that in many cases 
the turnover occurs between ~5-7 keV. They demonstrate it could 
be due to a combination of the iron K a emission line and ab- 
sorption edge (6.4 and 7.1 keV respectively) in a relativistically 
blurred reflection spectrum from the inner accretion disc of a spin- 
ning black hole. Here 'r eflection' refers to the b ackscattering and 
fluorescence of X-rays l lGeorge&FabiarJll99ll) . In this case the 
excess soft emission often assumed to come directly from the ac- 



cretion disc is also due to blurred reflection, so the disc is not ac- 
tually observed, hence the inner temperature cannot be measured 
and the mass must be estimated by some other method. This ex- 
planation for the soft emission is conceptually similar to that pro- 
posed bv ,Goncalves & So ria (2006j). Other explanations have also 
been proposed, often based on emission from the aforementioned 
slim discs in which advection of radiation fro m the inner disc re- 
duces the observed luminosity from this region l lAbramowicz et all 
1988), although these have not been as successful in reproduc ing 
the observed spectra (see e.g. attempts in iGladstone et"ai]l2009l) . 

It is important to determine the origins of this curvature, as it 
is one of the properties of ULXs that distinguishes them from their 
less luminous XRB cousins. Here we present an XMM-Newton ob- 
servation of a new ULX candidate located in NGC 4517 (hereafter 
NGC 4517 ULXl), which was brought to light during work on a 
new catalogue of ULX candidates in nearby galaxies (Walton et 
al. 201 1, submitted). Both reflection and Comptonisation interpre- 
tations are applied during the spectral analysis and the models di- 
rectly compared. The paper is structured as follows: section |2] de- 
tails the data reduction, section [3] describes the spectral analysis, 
section|4]discusses and compares the two models, focusing on how 
they may be distinguished, and finally section |5] presents our con- 
clusions. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Host Galaxy 

NGC 4517 (also known as NGC 4437) is an Sc type spiral galaxy 
located at RA = 12''32'"45.6^ Dec = +00°06'54.0" with a 
distance of ~18 Mpc (assuming Ho — 73 km s^^ Mpc^^, 
^matter = 0.27, r^vacuum = 0.73; the redshift of NGC 4517 is 
z — 0.003764)^, and has an apparent edge-on orientation. Al- 
though it has received relatively little observational attention, its 
orientation has lead to NGC 4517 being a popular galaxy to in- 
clu de in studies of Glob ular Clusters (e.g. iGoudfrooii et alj[2003l 
and IChandar et alj|2004l) . pot entially important e nvironments for 
the formation of IMBHs (see iNovola et alj|2008l . and references 
therein ) and hence possibly ULXs. Ba s ed on the diagnostic sys- 
tem of lKewlev eraU | |2006|) . ISeth et alj ( |2008|) find nuclear emis- 
sion line ratios consistent with star formation being a significant 
source of energy within NGC 45 17. Star for ming regions are als o 
important environments for ULX formation jSwartz et al.ll2009h . 
However, they also find evidence of AGN activity, and do not go as 
far as to provide a quantitative estimate for the star formation rate. 

The only targeted X-ray observation of NGC 45 17 was 
taken in December 2004 when it was observed with XMM- 
Newton jjansen et alj 1200 ih for ~114 ks, during which a num- 
ber of X-ray point sources were detected wit hin its D25 isophote, 
as de fined in the RC3 galaxy catalogue |de Vaucouleurs et alj 
Il99lh . including a point source at RA = 12''32'"45.0^ Dec 
= +00°06'55.0", co-incident with the expected location of the 
galaxy nucleus. However, the brightest of these sources, with 
an observed (absorbed) X-ray luminosity of Lo. 2-12.0 ~ 2 x 
10*° erg s~^ (assuming isotropic emission and that the association 
with NGC 45 17 is correct), is NGC 45 17 ULXl (lAU ID: 2XMM 
J123242. 7-1-000654), offset from the nucleus by 43 arcsec (see Fig. 

^ The basic information for NGC 4517 presented here has 
been obtained from the NASA Extragalactic Database (NED); 
|http://nedwww.ipac. caltech.edu/ 1 



© 2009 RAS, MNRAS QOQ.fTlflll 



Comparing Spectral Models for ULXs with NGC4517 ULXl 3 




Table 1. A comparison of the single and broken powerlaw models, modified 
by neutral absorption with column density Nh ~ 8 x 10^^ atom cm~^, 
applied to the 2-10 kcV data for NGC4517 ULXl (see text). 



8x10* 



Figure 2. The £P/C-pn lightcurve of NGC 45 17 ULXl with 1000 s resolu- 
tion, obtained during the XMM-Newton observation of NGC 4517 taken in 
December 2004. 

[B, with position RA = 12''32"'43.0", Dec = +00°06'55.0". It is 
this source that is the subject of this work. In contrast, the X-ray 
detection coincident with the galactic centre is observed to have a 
luminosity of only Lo. 2-12.0 ~ 10^^ ergs-\ The X -ray detec- 
tions of the ULX candidate and the nuclear source have position 
errors of 1.1" and 1.8" respectively (99 per ce nt confidence level ), 
as quoted in the 2XMM Serendipitous Survey dWatson et alj2009l) . 
while the NED position for the centre of NGC 4517 should be ac- 
curate to within 3.8", so we may be confident NGC 4517 ULXl is 
not associated with the nucleus of its host galaxy. 



2.2 Data Reduction 

The data reduction was carried out with the Science Analysis Sys- 
tem (SASv8.0. 1) largely according to the standard prescription 
provided in the online guide^. The observation data files were 
processed using EPCHAIN a nd EMCHAIN to p roduce calibrated 
event lists for the E PlC-pn jStruder et al1l200ll) and EPIC-MOS 
jTumer et al.l l20oi]) CCDs respectively. Lightcurves and spectra 
were produced for the 0.3-10.0 keV energy range selecting only 
single and double events (single to quadruple events) for EPIC-pn 
(EPIC-MOS) using XMMSELECT, and periods of high background 
were treated according to the method outlined by iPiconcelli et al.l 
( I2OO4I) with the signal-to-noise ratio maximised for the full con- 
sidered energy band. A circular region of radius 29"/6.7 pix. cen- 
tred on NGC 45 17 ULXl was defined in order to include as many 
source counts as possible without contamination from the nearby 
nuclear source or including CCD chip edges. The observation was 
targeted at the galaxy position, so NGC 45 1 7 ULX 1 was offset from 
the nominal pointing by ~45" for all detectors. A second, larger 
circular region of radius 99"/22.7 pix. was chosen in an area of 
the same CCD free of other sources to sample the background. We 
also extracted the spectrum of the source coincident with the nu- 
cleus, using a circular region of radius 16"/3.7 pix. so as not to 
overlap the region adopted for the ULX. The redistribution ma- 
trices and auxiliary response files were generated with RMFGEN 
and ARFGEN. Lightcurves were corrected for the background count 
rate using EPICLCCORR"'. After performing the data reduction sep- 
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arately for each of the MOS CCDs, the spectra were combined 
using the FTOOL** ADDSPEC to improve the statistics at the high- 
est and lowest energies. ADDSPEC combines spectra in a response 
weighted manner, to account for any differences there may be in 
e.g. the extraction regions etc., and also automatically combines the 
instrumental responses and background spectra associated with the 
source spectra. Finally, spectra were re-binned using GRPPHA to 
have a minimum of 25 counts in each energy bin, so that the prob- 
ability distribution of counts within each bin can be considered 
Gaussian, and hence the use of the statistic is appropriate when 
performing spectral fits. 

This procedure yielded good quality spectra for NGC 45 17 
ULXl, comparable with some of the best ULX spectra cur- 
rently obtained. In total, '--^ 15500 counts were recorded with the 
EPIC cameras, with ~8500 of these recorded by EPIC-pn, so the 
data m eet the quality criterion adopted by both [Gladstone et al.l 
( l2009h and lStobbart et alj ( |2006|) . The total good observation times 
contributing to the reduced spectra are ~70ks for EPIC-pn and 
~85 ks for each EPIC-MOS CCD. The EPIC-pn lightcurve ob- 
tained for the first 80 ks of the observation (of which ^^70 ks has 
a low enough background count to be considered good time and 
contribute to the spectrum) is shown in Fig. (2] during which the 
source remains fairly constant. After this time, the observation is 
dominated by very strong background flarin g events. We attemp t 
to calculate the fractional excess variability jEdelson et al.ll20()2h . 
but find that the lightcurve is completely dominated by statistical 
fluctuations and does not display any observable intrinsi c variabil- 
ity, sim ilar to a number of other ULXs highlighted in iHeil et al] 
( l2009h . 



3 SPECTRAL ANALYSIS 

Here we present our spectral analysis of NGC 45 17 ULXl. 
Throughout this work, spectral modelling is performed with 
XSPEC vl2.5.1n ( Arnaud 1996), and parameters are obtained by 
modelling the EPIC-pn and combined £P/C-MOS spectra simul- 
taneously. All parameters are tied between the two spectra with 
the exception of a constant multiplicative component, required to 
be 1.0 for EPIC-pn but allowed to vary for £P/C-MOS to account 
for possible cross-calibration uncertainties. In all cases presented in 
the following sections this parameter is within 5 per cent of unity. 
All quoted uncertainties are the 90 per cent confidence limits for a 
single parameter of interest, unless stated otherwise. The Galactic 
absorption c olumn in the direction of NGC 4517 is ~L88 x 10^" 
atomcm"^ jKalberla et al.i,2005i) . 



^ http://xmm.esac.esa.int/ 

^ Note that EPICLCCORR only accounts for the instrumental good time 
interval (GTI) generated during the production of the events file, it does not 
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account for additional GTI files based on e.g. periods of high background, 
as used here. 
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Figure 3. The EPIC-pa (black) and combined EPIC-MOS (red) 2-lOkcV spectra of NGC4517 ULXl, modelled with absorbed single and broken-powerlaw 
models (left and right panels respectively, see text); data/model ratio plots are also shown for each. The single powerlaw model clearly predicts an excess of 
counts over those observed at high energies (above ~6 kcV). The data have been rebinned for display purposes only. 



3.1 High Energy Turnover 

As previously stated, high energy curvatur e is frequently seen in 
the best quality X- ray spectra of ULXs ( [Gladstone et alj [20091 ; 
[Stobbart et al.[[20()6l) . Given that the data presented here is of simi- 
lar quality, it is natural perhaps to expect that such a turnover should 
also be observed in NGC 45 1 7 ULX 1 . Accordingly, we investigated 
the high energy data to determine whether a turnover was present 
by comparing single and broken powerlaw models to the 2- 10 keV 
spectra, following a similar prescription to [ Stobbart et al.[ l [2006l) . 
However, inspection of the full 0.3-10.0 keV data for NGC 4517 
ULXl suggests there is a significant amount of neutral material 
obscuring the source, with a column density of A^h — 8 x 10^^ 
atom cm^^ (far in excess of the Galactic column). This is most 
likely due to the edge-on orientation of the host galaxy, and unfor- 
tunately will add to the difficulty in reliably determining the pres- 
ence and origins of any excess emission at soft energies, and hence 
to the difficulty in distinguishing between spectral models with the 
current data. Modelling the (admittedly very poor) spectrum of the 
source associated with the nucleus with a simple absorbed power- 
law also implies significant absorption, with a column density of 
Nh^2x lO^^atomcm" ^. 

While [Stobbart et al.[ feOOd) compare unabsorbed models, the 
column required here will have some effect on the spectrum above 
2 keV, so we include absorption with a column density A^h — 
8 X 10^^ atom cm^^ in our comparison. We find that the broken 
powerlaw model is significantly preferred, providing an improve- 
ment of Ax^= 77 for 2 additional degrees of freedom (d.o.f.), the 
probability of chance improvement given by an F-test is negligible, 
at 2 X 10"". The comparison is shown in Table [T] and Fig. [3] in 
which it is clear there are fewer counts observed at high energies 
than predicted with the single powerlaw model. This implies that 
the high energy spectrum of NGC 45 17 ULXl does indeed display 
curvature, similar to that seen in other high quality ULX data. In 
addition, the break energy of 5.5 keV obtained is similar to that 
commonly seen in the other high quality ULX spectra, where the 
break energy is often close to ~6 keV. 

Having demonstrated the presence of intrinsic high energy 
curvature, in the following sections we apply the physically mo- 
tivated disc reflection and Comptonisation interpretations proposed 
to explain this feature. We do not preclude that other interpretations 
for this source are possible, but we limit ourselves to these two as 
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Figure 4. The disc reflection model applied to the spectrum of NGC 45 17 
ULX 1 and the relative contributions of the components included; the solid 
black line is the total model, the red shows the powerlaw continuum and the 
blue shows the blurred reflection component. Only the EPIC-pa data are 
shown for clarity. 

a demonstrative process in order to investigate how they may be 
distinguished in general. 

3.2 Disc Reflection Interpretation 

We begin by considering the reflection model presented in 
[Caballero-Garcia & FabiaiJfeOlOh . in which the observed turnover 
is due to the combined presence of a relativistically blurred iron 
emission line and absorption edge. These features arise from X-ray 
reflection within the inner regions of the accretion disc after irra- 
diation of this material by the intrinsic po werlaw continuum. T his 
model makes use of the REFLIONX grid | [Ross & Fabian[[2005l) . a 
self consistent reflection model which intrinsically includes iron 
K-shell absorption and emission. Its key parameters are the iron 
abundance of the reflecting medium (relative to the solar iron abun- 
dance), Apc, the photon index of the assumed ionising powerlaw 
continuum, F, and the ionisation parameter of the surface of the 
reflecting medium, ^ — L/nR^ , where L is the incident luminos- 
ity able to ionise hydrogen, n is the number density of hydrogen 
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Figure 5. Confidence contours of the inclination paired with the other blurring parameters: i?in, dinner, Rhr and goutor (left to right). The contours shown 
are the 90, 95 and 99 per cent confidence levels for two parameters of interest (black, red and green respectively). In each case the best fit parameter values 
given in Table |2] are marked with a + symbol. There do not seem to be large parameter degeneracies present in the reflection interpretation presented for 
NGC4517 ULXl at even the 95 per cent level. 



and R is the distance from the ionising source. REFLIONX does 
not include the abundances of any elements other than iron as free 
parameters, the m odel is instead calculated assum ing their solar 
values as given in iMorrison & McCammoij ( 1 19831) . The relativis- 
tic blurring is applied using the convolution model KDBLUR which 
accounts for the extreme gravitational blurring expected in the in- 
nermost r egions aroun d a Kerr black hole using calculations per- 
formed bv lLaoj 1 IT991I) . Here the key parameters are the inner and 
outer radii of the accretion disc, Rin and -Rout respectively, its in- 
clination with respect to the observer, i, and its emissivity index, 
q, which describes the assumed powerlaw emissivity profile of the 
disc, of form e(r) cc r^'. Higher values of q imply preferential il- 
lumination of the inner regions of the disc, and hence a more com- 
pact, centrally located corona (we caution the reader that although 
there are differences in the Compton scattering medium between 
the reflection and Comptonisation interpretations, in terms of the 
physical condition and possibly also the origin of the matter, as 
will be discussed later, we refer to it as the corona in both cases). 

Initially we applied the mode l following the same method as 
ICaballero-Garcfa & Fabiai] ( 1201 Ol) including both Galactic and in- 
trinsic neutral absorption using PHABS and ZPHABS respectively. 
A reasonable fit with xS(d.o.f) = 1.2(468) is obtained with sim- 
ilar parameters to those presented in the initial paper, i.e. signifi- 
cant blurring of the reflected emission, with the accretion disc ex- 
tending to the last stable orbit of a maximally rotating black hole. 
The disc is required to have a highly super-solar iron abundance, 
and the ionising continuum is extremely steep in this case. How- 
ever, the turnover was still not adequately modelled. To resolve 
this, we make a minor modification to the model (without changing 
the physical components included) and replace KDBLUR with KD- 
BLUR2, improving the fit to x^(d.o.f) = 1.1(466) (Ax^of 40 for 2 
extra d.o.f.). KDBLUR2 applies the relativistic effects due to strong 
gravity in the same way as KDBLUR, but allows the emissivity pro- 
file of the disc e{r) to have a broken powerlaw form, breaking at 
-Rbr- The parameters obtained with this modification are given in 
Table|2l and the relative contributions of the components are shown 
in Fig.|4] We still find the emissivity profile to be strongly centrally 
peaked. During the application of the reflection model, the outer 
radius of the disc was set to 400 Rg, the maximum allowed by the 
KDBLUR and KDBLUR2 models, as it is very poorly constrained 
due to the centrally peaked emissivity index. The only constraint 
obtained if this parameter is allowed to vary is that the outer radius 
of the disc must be greater than 10 Rq. 

The disc reflection interpretation constructed is a complex, 
multi-parameter model. We therefore investigated the confidence 
contours for various combinations of parameters to search for pa- 
rameter degeneracies, focussing on those of the KDBLUR2 compo- 



Table 2. Parameters obtained with the reflection interpretation for the spec- 
trum of NGC 4517 ULXl; parameters marked with * have not been allowed 
to vary. 



Component 


Parameter 


Value 


PHABS 




0.188* 


ZPHABS 








z 


0.003764* 


POWERLAW 


r 


> 3.18 


KDBLUR2 


ifin" 


-1 9^9+0.065 






400* 




i 






•Jinncr 


> 9.18 




•Joutcr 


4 ^6+0'^^ 




Rbr'' 




REFLIONX 




1 3+^''^ 
^■■'-0.3 






> 8.2 


X^(d.o.f) 




1.1(466) 


-I'D. 3-10'' 




24.611;? 



" Column densities are given in lO^i atom cm ^ 
Radii are given in units of gravitational radii, Rq = GMgfi / c? 
lonisation parameter, given in erg cm s~i 
iron abundance, quoted relative to the solar value 
The absorption corrected 0.3-10.0 koV luminosity in 10^" erg s^l 

nent. Degeneracies in the blurring parameters have been found to 
be a significant issue in a number of cases d uring the applicatio n 
of reflection models to other sources (see e.g. lNardinietal]|2010l) . 
Fig. [5] shows the confidence contours of the disc inclination paired 
with each of the remaining KDBLUR2 parameters (excluding i?out, 
which was not free to vary). It is clear from these panels that strong 
parameter degeneracies are not an issue for the reflection interpre- 
tation of this particular source, even at the 95 per cent confidence 
level for two parameters of interest. 

3.3 Comptonisation Interpretation 

The second model applied is that of [Gladstone et all ( |2009|) . which 
proposes a new accretion state for X-ray binaries, dubbed the 'Ul- 
traluminous State' . In this scenario, the turnover is due to Compton- 
isation of disc photons by an optically thick corona. To model this. 
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Figure 6. The Comptonisation model applied to the spectrum of NGC4517 ULXl, and the relative contributions of the components; the outer (observable) 
disc is shown in red and the Comptonised emission in blue, while the overall model is in black. Again only the EPIC-pn data are shown for clarity. The left 
panel shows the model fit to the photon spectrum, while the right panel shows the proposed model (absorption corrected) and its components in an EFe plot, 
highlighting the effects of the corona on the observed disc emission; the additional emission from the disc that would be observed if the corona and its effects 
were removed is shown in green. 



we make use of the DKBBFTH code jPone & Kubot j|2006h . which 
accounts for the e nergetic couphng between the disc and the corona 
using the model of lSvensscm &Zdziarskil l ll994l) . i.e. it assumes the 
corona extends homogeneously (constant temperature and optical 
depth) over the disc out to some transition radius, Rc. The corona 
takes some fraction, /, of the gravitational energy available within 
Rc, with the remaining fraction (1 — /) powering the inner disc 
emission. As a simplifying assumption the model also takes / to be 
constant with radius. Other authors have prop osed alternative sce - 
narios in which / may vary wit h radius (e.g. Ijaniuk et aLlboOOl) . 
but as shown in iDone & Kubotal l l2006h . keeping / constant is an 
adequate approximation for even the best quality XRB data, and so 
will be sufficient for our needs. The key parameters of DKBBFTH 
are the temperature of the Comptonising electrons, Tc, the spectral 
index of the corona, F, the transition radius (or the outer radius of 
the corona), and the inner temperature of the accretion disc, Tin. 
Note that this temperature is what would be observed if the corona 
and its effects on the disc were completely removed. The fraction of 
the power in the corona / is calculated self consistently via itera- 
tion from the coronal spectral parameters. The inner radius of the 
disc, Tiin, may be calculated from the normalisation of the model 
if the distance to the source and the inclination of the (inner) disc 
are known, in the same manner as for the DISKBB model, and the 
optical depth, r, can be calculated from F and Te (assuming a slab 
geometry). 

We applied this model to the data, allowing for both Galac- 
tic and intrinsic neutral absorption in the same way as for the re- 
flection interpretation, using PHABS and ZPHABS again, and ob- 
tainedagood fit with xt{d.o.i) = 1.1(471). Model parameters, and 
key quantities derived from them, are given in Table [3] and the fit 
to the data are shown in Fig. |6] As with the reflection model, we 
fi nd similar parameters h ere to those obtained in the original work 
of [Gladstone et al.l ( l2009h . The electron temperature of the Comp- 
tonising corona is low (fcTc ~ 1.4 keV), and it is optically thick 
(r ~ 9). In addition, the temperature of the innermost disc inferred 
from the model (fcTin ^ 0.9 keV, in the limit of no corona) is akin 
with that seen in stellar mass black hole binaries in accretion states 
with high Eddington ratios, e.g. the high/soft state or the steep pow- 



Table 3. Parameters obtained, and quantities derived from them, with the 
Comptonisation interpretation for the spectrum of NGC 4517 ULXl; pa- 
rameters marked with * have not been allowed to vary. 



Component 


Parameter 


Value 


PHABS 


Nn" 


0.188* 


ZPHABS 


Nh" 


4 2+0-4 
^■^-0.3 




z 


0.003764* 


DKBBFTH 


r 


1 7fi+"'09 

J^- "'-o.oe 






1 4+0.2 






9+°-^ 

U.M_g 2 




RinV COS i 


ogn+TOO 

•^^"-250 




Rc/RiJ 








9 04+0'09 


xUd.o.{) 




1.1(471) 


io.3-l(/ 




3.53 ±0.05 



° Column densities are given in 10^^ atom cm~^ 

^ Temperatures are given in kcV 
Inner radius of the disc in km. in terms of the unknown incUnation angle 
Outer radius of the corona, relative to the inner radius of the disc 
T is not a model parameter, but is calculated from F and Te 

f The absolution connected 0.3-10.0 keV luminosity in 10'"' erg s~^ 



erlaw state (using the nomenclature of iRemillard & McClintoc^ 
l2006h . The inclination angle of the inner disc, i, is not known in 
this case, so the inner radius of the disc is quoted in terms of its 
functional dependence on this angle. 



4 DISCUSSION AND COMPARISON 

We find that both disc reflection and Comptonisation can model the 
high quahty 0.3-10.0 keV EPIC spectra of NGC 4517 ULXl wefl 
when a broad parameter range is considered, and in fact give statis- 
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tically equivalent fits to the data. This appears to be the case for the 
majority of ULXs wit h spectral curvature, as all the sour ces mod- 
elled with reflection in ICaballero-Garcia & Fabiaiil ll20ld) are also 
well modelled with Comptonisation bv iOladst one et alj j2009l) . 
Currently these spectra represent the best available ULX data, so 
other means are necessary to distinguish between these two inter- 
pretations. We begin by considering the physical implications of 
our analysis, as well as the results obtained for other ULXs. 

The lack of direct thermal emission from the accretion disc in- 
cluded in the reflection model may at first seem puzzling. For many 
of the black hole masses invoked in both the stellar mass and IMBH 
interpretations, the inner temperature of the disc is expected to oc- 
cur in the X-ray bandpass (~1.0 keV for a ~10 M0 stellar mass 
black hole, down to ~0.3 keV for a ~1000 M© IMBH). Given 
the strong disc reflection, one might expect to see a sizeable con- 
tribution from a direct thermal component. However, if magnetic 
extraction of energy from the di sc to the c orona (which itself might 
be m agnetic in nature, see e.g. [Merloni & Fabian 2002, Liu e t akl 
I2OO2I and references therein) is highly efficient, energy may be ex- 
tracted from the disc before it is able to thermalise. In this way, the 
direct thermal emission could be suppressed with respect to that 
of the corona, and hence also the reflected component. Such pro- 
cesses are not currently well understood and merit further attention. 
Were this to be the case, such highly efficient magnetic extraction 
would distinguish ULXs from both XRBs and AGNs, as it is well 
known that there is often evidence for direct thermal emission in 
both cases (albeit in the UV for AGNs), especially at high accre- 
tion rates (as are anticipated to be present in ULXs). Hence the 
efficiency of magnetic extraction is not expected to be particularly 
high in XRBs and AGNs. 

Of course, the direct thermal disc emission must be distin- 
guished from any of the reflected emission thermalised in the disc 
during the reflection process. The temperature at which this feature 
appears, for a given ionisation state, depends on the rate of free- 
free absorption and hence on the density of the illuminated surface 
of the disc. It is incorporated into the REFLIONX model, however, 
this model is calculated assuming a hydrogen number density of 
riH = 10^^ cm^'^, for which the thermalised component of the re- 
flected emission occurs at <0.1 keV; if the assumed density is not 
appropriate for the discs around ULXs then this emission will not 
be accounted for correctly. For NGC 45 17 ULXl this might not be a 
significant issue as the large absorption column prevents a detailed 
study of any additional emission at soft X-ray energies. For sources 
with less absorption it will be necessary to include this effect with 
reflection models calculated for a higher density (e.g. REFHIDEN; 
IRoss & Fab ian 2007,). This will be explored in future work. 

The emissivity profile obtained is strongly centrally peaked, 
which suggests that gravitational light bending is important. For 
this to be the case, the corona must be compact and located close 
enough to the central black hole that the gravitational potential 
focuses additional flux from the corona away from the observer 
and onto the the accretion disc. This allows the reflected compo- 
nent to ap pear more luminous than the intrinsic component; see 
iMiniutti & Fabian ( 2004) for a more detailed discussion of this 
phenomenon. In their work, it is shown that gravitational light 
bending preferentially illuminates the inner disc, and hence in such 
a regime the radial emissivity profile of the inner disc is not best de- 
scribed by a single uniform powerlaw, so the significant improve- 
ment obtained by allowing for a broken powerlaw emissivity profile 
is not unexpected. Given that centrally peaked emissivity profiles 
appear to occur fairly frequently when applying a disc reflection 
interpretation to ULX spectra, and that some sources require reflec- 



tion dominated solutions, allowing for a more complex emissivity 
profile may lead to similar improvements in the modelling of other 
sources. 

The proximity of the X-ray emitting corona to the black hole 
required to provide the obtained emissivity profile for NGC 45 17 
ULXl suggests that, for a corona on-axis, the majority of the flux 
ends up crossi ng the event horizon ra ther than illuminating the disc. 
The model of iNiedzwiecki & Zvckil ((2008.) requires the corona to 
be within ~1.5 Rq, at which point, if on-axis, ^--^90 per cent of the 
coronal flux may not illuminate the disc, so the true intrinsic lumi- 
nosity would be substantially larger than that observed. In addition, 
any small change in the height of the corona above the black hole 
would lead to large changes in observed luminosity, which may be 
in contrast with the e xtremely stable ULX lightcurves seen here and 
in Eeiret^ |200i). A more appropriate g eometry may be simi- 
lar to that proposed for MCG -6-30-15 bv INiedzwiecki & Zvckil 
( l2008l) . where the corona is located just above the inner disc. In 
this scenario the illumination from the corona is Doppler beamed 
along the disc, so a greater proportion illuminates the disc rather 
than being lost over the event horizon. 

At a first glance, the high iron abundance obtained for the ac- 
cretion disc her e, and in general when modelling ULX spectra with 
disc reflection jCaballero-Garcia & Fabian 201^, does not seem 
consistent with the observati on that ULXs are often located in low 
metallicity regions (see e.g. ISoria et al]|2005l iMapelli et alj|2009l 
iMapelli et alKOlOl) . nor the theoretical considerations that suggest 
the A/bh >10 Mq black holes ULXs are speculated to host are 
easier to form in low metallicity regions, both due to the formation 
of larger p rogenitor stars, and reduced mass loss rates from stel - 
lar winds jMadau & Reesll200ll ; lBromnjl2004l : [vink et alJl26oih . 
However, if the black holes in ULXs are accreting from debris discs 
formed from the metal rich inner layers of the stars that originally 
collapsed to form them in the first place, it may be possible for the 
accretion disc to have a very high metallicity even though the sur- 
rounding regions do not. Such discs could be produced by a method 
similar to that proposed by Li (2003), where the inner regions of 
the collapsing/exploding star are not properly ejected, and fall back 
towards the collapsed core to form the accretion disc. In such a sce- 
nario the black hole need not be in a mass transfer binary system 
to appear bright in X-rays, and we note that a number of ULX op- 
tical counterparts appear more consistent with X-ray irradiation of 
the outer regions of an acc retion disc than with a stellar companion 
(e.g. iKaaret & FenH2009l) . although it may just be that the emis- 
sion from the disc dominates that of any companion star even in 
the optical bandpass, as is the case for LMXBs in outburst. Debris 
discs might be a common scenario if ULXs are associated with su- 
pernova remnants (SNR), an explanation initially proposed for the 
ionised nebulae man y ULXs are embedded wit hin, e.g. IC342X-1 
([Roberts et al.ll2003l) and MF 16 in NGC 6946 ( [Roberts & Colberj 
I2OO3I) . 

However, there is growing evidence that these nebulae are not 
SNR, but a combination of X-ray photoionised and shock-ionised 
material due to irradiation and outflowing winds from the central X- 
ray source (see e.g. [Pakull & Grise lOOS). From these nebulae the 
ages of ULXs can be estimated as ~10''"^ yr dPakull et al1l2006l; 
IPakull & Gri"sil2008l) . which can be used to estimate the mass of 
high metallicity material required in the debris disc scenario. In 
order to radiate at I/x > 10'^^ erg s~^ the black hole must be ac- 
creting at AI > 10"'' Mq yr~^, so in order for ULXs to have 
observable debris discs, these must have formed from at least --^0.1 
Mq of high metallicity material. Of course, this makes the simplis- 
tic assumption that the mass accretion rate is constant with time. 
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iMineshige et alj l ll997h studied post supernova accretion via fall- 
back of matter in detail and found that the accretion rate should 
evolve with time, as M oc with 1 < a < 1.5 (the exact value 
of a depends on the conditions of the material as it begins to fall 
back). In this case the required accretion rate of M > lO"'^ M© 
yr~^ would represent the 'current' accretion rate, and the initial 
mass of high metallicity material available may have had to have 
been much greater than 0. 1 M© even for a ULX currently radiat- 
ing at 10^^ erg s^^. In addition, disc winds may have been present 
throughout, or at stages during the lifetime of the ULX. The initial 
disc masses required are therefore at least an order of magnitude 
larger than the lower limits for the current disc mass estimates for 
th e ULXs CXO J033831.8-352604 and XMMU122939.7-I-075333 
in |Porteil ( l2010l) . 

It may be difficult to form such large discs through fallback. 
Supernovae in low metallicity regions may not eject large amounts 
of material jZampieri & Robertal2009l) . hence massive, rapidly ro- 
tating progenitor stars may be required. In addition, a number of 
optical ULX counterparts have been identifie d as probable high 
mass or evolved stellar companions, (see e.g. Kuntz et al. 
ISoria et all I2OO5I. iLiu et all l2007l [Roberts et all l2008l but espe- 
cially Patruno & ZamDierill2008h , from which mass is expected to 
be transferred to the black hole. The binary companion and the 
mass transfer could disrupt and dilute the metallicity of such a 
large disc. The possible presence of outflows is also interesting. 
The population of optical nebulae associated with ULXs appear to 
sho w a combination of ph otoionised and shock-ionised material, 
and IPakull & Grisi ( l2008h argue the latter cases may be inflated 
by winds from the central source, although the possibility that they 
are remnants of the explosive event in which the ULX was formed 
is not ruled out. If these are due to winds and ULXs are accreting 
from high metallicity fallback material, some metal enrichment of 
the regions surrounding the ULX might be expected, although ex- 
actly how much effect on the elemental abundances these outflows 
should have will depend on the amount of material ejected in com- 
parison to that in the ISM. It seems that although debris discs can 
not be ruled out for a fraction of the population, they may not be a 
viable explanation for all ULXs. 

An alternative explanation is that the iron abundances obtained 
may not be good estimates of the true abundances of the systems 
in question, but systematically modified by processes that are not 
accounted for in the modelling. For example, this would be the 
case if there are actually multiple reflections of the intrinsic con- 
tinuum. The observed emissivity profile implies that strong gravi- 
tational light bending is an important process in NGC4517 ULXl, 
and ULXs in general, if reflection is the correct interpretation. This 
causes the initial reflection of the intrinsic powerlaw component to 
occur primarily within the central regions of the accretion disc, so 
the majority of the reflected emission also originates within the re- 
gion in which strong light bending is present. Some fraction of this 
reflected emission could therefore be bent back towards the disc 
and reflected a second time. The initially reflected emission already 
has the iron emission features imprinted onto it and the secondary 
reflection should serve to enhance these further, as demonstrated by 
IRoss et alj ( l2002h . hence the iron abundance obtained from mod- 
elling the emission with a single reflection process could be artifi- 
cially inflated. If this is the case, it may not be necessary to invoke 
the high metallicity debris discs discussed previously. Whatever the 
correct explanation for the super-solar iron abundances obtained, 
it might be natural to expect that other heavy elements might be 
or appear over-abundant for the same reason. Unfortunately model 
limitations currently prevent us from exploring this possibility. 



The Comptonisation interpretation requires that the coro- 
nal electron temperature of NGC4517 ULXl is rather low, at 
fcTc ~ 1 keV. Indeed, low electron temperatures (fcTe ~ 2 keV) 
are a c ommon feature in the application of this model to ULX 
spectra ( [Gladstone etai]|2009l) . as they are an unavoidable con- 
sequence of modelling the observed curvature as the high en- 
ergy cut-off of Compton up-scattering by thermal electrons. Typ- 
ical electron temperatures for XRBs are usually much higher, 
fcTc ~ 50 — 100 keV for binaries in the low/hard state and 
~30-50 keV in the high/soft state (although here there can also 
be a significa nt contribution from a non- t hermal electron distribu - 
tion); see e .fj. 'Malza c & Belmonll | |2009|) . iGierlinski et"ai] ( Il999t) . 
ISunvaev & Titarchuk l ll980h amongst others. The coronal opti- 
cal depth is also quite discrepant from typical values for XRBs: 
r~ 0.1 — 0.3 in the high/soft state and ~1 in the low/hard state. 
Here we find r ~ 9. However, ULX spectra are often likened 
most to the steep powerlaw state. Sources in this state often dis- 
play very high Eddington Ratios (L/Le ~1). Again, Compton up- 
scattering by thermal electrons is observed to be important, with 
fcTc - 10 - 30 keV llWilson & Don"3l200ll;lKubota & Makishimal 



I2OO4I : lKubota& Donell2004h . These temperatures are often still 
highly discrepant from those obtained from ULX spectra, although 
the highest luminos ity observations of the XRB GRS 1915-1- 105 are 
notable exceptions jMiddleton et alj2006l : IUeda et al.ll2009l) . 

However, in the work of [Gladstone et"ai] ( |2009|) . it is sug- 
gested that most ULXs may represent a new accretion state en- 
tirely, dubbed the 'Ultraluminous State', which sources enter if 
and when they are able to radiate at super-Eddington luminosi- 
ties (I//I/E > 1). As sources approach, and potentially exceed the 
Eddington limit, there is much speculation that r adiatively driven 
disc w inds become increasingly important (see e.g. lKing & PoundU 
l2003h . A case in point, potentially relevant for ULXs, is the bi- 
nary system SS433 which is expected to have a very high mass 
transfer rate and from which large outflows are observed, in addi - 
tion to relativistic jets /Beg elman et alj200d ; |Poutanen et al.l2007h . 
[Gladstone et al. (2009) argue that such winds may be the source 
of the increased coronal opacity obtained from ULX spectra, as 
further mass is injected into the corona by the outflow. The elec- 
tron density of the corona is then increased, and there is therefore 
less energy available per electron, so the average temperature of the 
electrons decreases. These outflows could be responsible for the in - 
flation of the observed ULX optical nebulae l lPakull & Grisd2008h . 

In brief, it seems that on purely physical grounds, neither in- 
terpretation can be excluded from being the underlying process for 
at least some fraction of the ULX population, so we must seek 
other, observational means to distinguish between them. Ideally, 
this would be easily applied to a sample of ULXs large enough to 
be considered representative of the population as a whole. Unfor- 
tunately, the models do not make unique predictions for the emis- 
sion at non X-ray wavelengths, so multiwavelength studies are un- 
likely to resolve this particular issue. In Fig. |7] we compare the 
two models fit to the spectrum of NGC4517 ULXl, but extended 
beyond the XMM-Newton bandpass up to 100 keV. A few minor 
details aside, the models agree very well below ~10 keV. As we 
have demonstrated an observation of ~100ks is insufficient to dif- 
ferentiate between them for NGC 45 17 ULXl, and the wo rk of 
ICaballero-Garci'a & Fabiar] ( 1201 Ol) and lGladstone et al.l ( |2009|) sug- 
gests both models agree well at these energies for similar length ob- 
servations of even brighter ULXs. Attempting to determine which, 
if any, of these models is correct using only CCD data would re- 
quire even longer observations, which are not realistically obtain- 
able for a large number of sources. Some ULXs may be bright 
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Figure 7. A comparison of the Comptonisation and disc reflection mod- 
els for ULX spectra applied to NGC4517 ULXl; both have been extended 
beyond the XMM-Newton bandpass up to 100 koV. Although the models 
agree well below lOkeV, there is a clear discrepancy of many orders of 
magnitude between them at higher energies, due to the presence of the 
Compton hump in reflection spectra. Inset: data simulated using the cur- 
rent Astro-H response files for a lOOks observation with both models. 
NGC4517 ULXl is clearly detected to higher energies with the reflection 
interpretation. 



enough that the reflection grating spectrometers aboard XMM- 
Newton could investigate the differences in the neutral absorption 
the two models predict with reasonable observation lengths. How- 
ever, these are very much in the minority, and obtaining an inde- 
pendent value for the column density would not confirm either in- 
terpretation directly, it would merely provide an an additional con- 
straint that each model must fulfil if they are to continue to provide 
acceptable representations of the data. 

Potentially, the most fruitful avenue utilising only data below 
10 keV is the combination of spectral and variability studies. In- 
vestigation of the ms variability spectrum can aid significantly in 
the deconvolution of separately varying spectral components and 
could be able to differentiate between the two models considered 
here. Both are two component solutions, but the energy depen- 
dence of the relative importance of the two components is different 
for the two models. This approach is being investigated by Mid- 
dleton et al. (in preparation). However, a reliable estimation of the 
rms variability spectrum requires good quality data, which often 
requires long observations. In addition, this me thod relies on the 
ULXs displaying observable variability, however [Heil et ah 
demonstrated that the brightness of many ULXs is remarkably sta- 
ble within single XMM-Newton observations up to ~ 100 ks in dura- 
tion; it seems that NGC4517 ULXl may also be counted amongst 
this number. As a further complication, there are also sources that 
exhibit flux variabilit y without any strong spectral variability (see 
iKaaret & Fen3l2009l ) This might be considered unusual for inter- 
pretations invoking two distinct emission components. However, 
these components are not independent of one another, in both mod- 
els one component represents reprocessed emission from the other. 
If the luminosity of the primary component varies, then the lumi- 
nosity of the secondary component must also vary in the same way. 



A brightening of the primary component might lead to some vari- 
ation in the observable parameters, for example the temperatures 
of the observable part of the disc and the Comptonising electrons 
in the Comptonisation interpretation or the ioni sation of the disc 
in the reflection interpretation, but as argued by IVierdavanti et al.l 
( I2OIOI) these effects are likely to be fairly subtle an d not easily de- 
tected with the hardness ratio analysis employed by lKaaret & Fend 
Unfortunately, combining spectral and variability studies 
to decompose emission components is often only possible when 
strong spectral variability is present. All of these requirements limit 
the number of sources for which such analysis is possible. 

Returning to Fig. |2l the comparison above ~10 keV shows 
the two models diverging. The Comptonisation model continues 
its downward curvature, as this is due to the high energy limit of 
the Compton up-scattering process of disc photons by cool elec- 
trons, while the reflection spectrum turns back up above 10 keV 
due to the Compton hump. This is a broad emission feature present 
in reflection spectra that peaks at ~20-50 keV due to the interplay 
within the reflecting medium between Compton down-scattering 
of high ener gy photons and photoel ectric absorption of low en- 
ergy photons jGeorge & Fabianll99lh . Features consistent with the 
Compton hump are fairly c ommonl y observed in the spectra of both 
BHBs and AGNs (see e.g. iReis et alj|201 0; Wa lton et al. 20ia) . The 
difference in the origins of the observed curvature between the two 
models leads to large differences in their predicted high energy 
spectra. Of course, the Comptonisation model assumes the elec- 
trons have a purely thermal distribution. It is likely there is an ad- 
ditional non-thermal electron population that will eventually arrest 
the downward curvature predicted in this Comptonisation model, 
but we still anticipate clear differences between the two models at 
high energies, as for any reasonable non-thermal temperature dis- 
tribution these electrons would merely limit the downward curva- 
ture and provide a high energy powerlaw tail, not turn the spectrum 
back up again. In addition, if ULXs are an extremely high accretion 
rate extension of the Steep Powerlaw State {i.e. the Ultraluminous 
State), the emission from the non-thermal electrons may not be sig- 
nificant until energies > 100 keV. 

The strength of the predicted Compton hump is obviously de- 
termined by the parameters of the reflection model, notably the in- 
clination of the disc. Although we find the key parameters to be 
well determined in this case, parameter degeneracies may exist in 
other cases. REFLIONX is an angle averaged model, best suited 
for use with inclination angles ~45°, so we investigate how the 
strength of the Compton hump may change for extreme inc lina- 
tions with the PEXRAV model jMagdziarz & Zdziarskilll995l) . At 
low inclinations (10°) its strength is actually enhanced by ~5-10 
per cent (relative to the 45° case), while at high inclinations (80°) 
its strength is reduced by ~25-30 per cent. This variation is highly 
unlikely to be able to reconcile the predictions of the two mod- 
els in any given case. Other factors may also modify its strength. 
For example, REFLIONX is calculated assuming the illuminating 
spectrum is a powerlaw with a high energy cut-off at 300 keV. 
Lowering the energy of this cut-off will also serve to reduce the 
predicted Compton hump to some extent. However, we stress that 
when both models provide acceptable representations of the data 
below 10 keV, the extra emission due to the presence of this high 
energy feature in the reflection interpretation should lead to large 
differences in the predictions of the two models above 10 keV for 
just about any parameter combination. Therefore, observations of 
ULXs at such energies should be an extremely important diagnostic 
of the underlying physical processes powering the emission from 
ULXs. 
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Table 4. Key information for the Hard X-ray Detector (HXD) and the Burst 
Alert Telescope (BAT) aboard the Suzaku and Swift satellites respectively. 



Instrument 


Type 


Energy Range 


Angular Resolution 






(keV) 


(FWHM, arcmin) 


HXD (PIN) 


Collimating 


12-70 


34 (FoV) 


BAT 


Coded Mask 


15-150 


17 (PSF) 



At the time of writing there are two facilities that are ca- 
pable of observing simultaneously in the ~0.5-10.0 keV range 
required to determine the key model parameters for each inter- 
pretation from the curvature at ~6 keV, and at the higher en- 
ergies required to d istinguish between them. These are Suzaku 
jMitsuda et al.ll2007l) . with the combination of the Hard X-ray De- 
tector (HXD, specifically the PIN detector; 'Takahashi et al."2007') 
and X -ray Imaging Spect rometers (XIS; |Kovam a et al. 2007), and 
Swift jGehrels et al]|2004l) with the combination of the Burst Alert 
Telesc ope (BAT; Barthelmv et alJlToOSi) and the X-ray Telescope 
(XRT; iBurrows et al.ll2005l) . Some key information on the hard X- 
ray instruments on board these two satellites is summarised in Table 
|4] Both have very modest spatial resolution (the HXD is a collimat- 
ing instrument, so its effective spatial resolution is the same as its 
field of view). For the vast majority of ULXs, and certainly for the 
specific case of NGC4517 ULXl we have been considering here, 
neither of these instruments would be able to resolve the ULX from 
the nucleus of its host galaxy, which in most galaxies (although not 
necessarily NGC 4517) will be the dominant source of hard X-rays. 
This has been a common problem with hard X-ray detectors to date, 
and as such very little is known about the spectra of ULXs above 

10 keV. 

The exception is M 82 X-1, for which iMivawaki et al.l ( |2009|) 
claim a Suzaku HXD detection up to '~-^20 keV. In this work, the 
3.2-20.0 keV spectrum of M 82 X-1 is found to require a curved 
continuum (below ^^3.2 keV the spectrum is dominated by diffuse, 
thermal emission). Some uncertainty remains over what propor- 
tion of the flux detected above 10 keV originates from M 82 X- 
1 and what proportion originates from the rest of the (sizeable) 
XRBAJLX population of M82, but the requirement of a curved 
continuum is not sensitive to this. This curvature is well modelled 
with Comptonisation, and similar parameters to the results from 
other ULXs are obtained (fcTc --^ 2.5 keV, r ^ 8). However, this 
was only possible due to the unique scenario M 82 X-1 is in. M 82 is 
a starburst galaxy with no discernible AGN activity, an d of its X-ray 
source population, M82 X-1 is usually the brig htest jZhang et al.l 
l2009l ; IZezas et al.ll200ll ; [Matsumoto etai]|200lh. In a ddition, dur- 



ing the observations utilised by Mivawaki et al. | |2009|) , it was pos 



sible to choose the pointing to ensure other nearby bright X-ray 
sources, e.g. the LLAGN in M 81, M 81 X-6 and HolmberglX X-1, 
were outside of the HXD field of view. 

It may be possible for the PIN and BAT detectors to contribute 
to the differentiation of other sources if they are also significantly 
brighter than any others that they can not spatially be resolved 
from. In addition, a total non-detection at a low enough flux sen- 
sitivity would rule out the reflection interpretation. However, the 
sensitivities of these instruments is rather limiting for the study 
of ULXs. We stress that these scenarios will both be extremely 
rare, so such an approach will not be possible for the vast major- 
ity of the ULX population. M 82 X-1 is the only currently known 
source these instruments have been able to study due to its appar- 
ently unique situation and high flux. In order to obtain spectral 
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Figure 8. Data simulated from the reflection interpretation presented for 
NGC 45 17 ULXl for the suite of instrumentation aboard the Astro-H satel- 
lite (black: soft X-ray imaging system; red: X-ray microcalorimeter; green: 
hard X-ray imaging system). The solid lines show the application of the 
Comptonisation model to this simulated data. Clearly there is an excess of 
counts over the model prediction at energies above ~ 10 kcV. 



information above 10 keV for a large sample of ULXs an imag- 
ing spectrometer with significantly higher spatial resolution than is 
currently available will be required. With the improvement in sen- 
sitivity and angular resolution (~60-90") offered by the Hard X- 
ray Im aging System proposed to fly on Astro-H^ jKokubun et al.l 
1200^ in 2014, combined with the Soft X-ray Imaging system and 
the Soft X-ray Spectroscopy System (microcalorimeter) also due 
to be on board, it should be possible to begin addressing the na- 
ture of the observed emission from ULXs with a series of relatively 
short, pointed observations. In the shorter term, such observations 
should also be possible wit h simultaneous observations from A';(5- 
TAR^ ( iHan-isonet al]|2005h . which is due for launch in 201 1, and 
either XMM-Newton or Suzaku. An excellent candidate with which 
to demonstrate the capabilities of these instruments to contribute 
to ULX studies is Holmberg IX X-1, which is bright, nearby, and 
has a very hard spectrum, but currently cannot be resolved from the 
nucleus of M81 above 10 keV. 

As an immediate demonstration of the potential of these two 
instruments for ULX science, we have obtained the pre-flight in- 
strument response files for A.stro-H and simulated the data that 
would be obtained for both the models obtained here with a 100 ks 
observation of NGC 45 17 ULXl (we only simulate for Astro-H be- 
cause the hard X-ray imaging responses are very similar for this 
mission and NuSTAR, and Astro-H will also carry soft X-ray spec- 
trographs). These simulations are shown in the inset in Fig. |7] 
NGC 45 17 ULXl would clearly be detected out to much higher 
energies with the reflection interpretation than with the Compton- 
isation interpretation. We apply the Comptonisation interpretation 
to the data simulated from the reflection model for the full suite 
of instrumentation aboard Astro-H, allowing its parameters to vary 
(although requiring that Comptonisation models the ^^6 keV cur- 
vature), and obtain a relatively poor fit, with xt (d.o.f ) = 1 .33(284). 
Even though the electron temperature of the corona increases to 
Tc ~ 3.6 keV and its optical depth decreases to t ~ 3.5, Fig. [8] 
shows there is a clear excess of counts at hard energies. To provide a 



http://astro-h.isas.jaxa.jp/si/index_e.html 
^' http://byl34w.bay 134. mail. live. com/default. aspx?wa=wsigninl.0 
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rough estimate of the number of ULXs for which Astro-H and M(5- 
TAR should be able to make such a distinction with observations of 
a similar or shorter duration to that considered above, there are 19 
sources in the ULX catalogue presented by Walton et al. (2011; 
submitted) with a greater flux than NGC 45 17 ULXl, including 5 
whose flux is more than an order of magnitude higher. We stress 
that this is actually a lower limit, as there are a number of bright 
ULXs that are not present in the catalogue, of which Holmberg IX 
X- 1 is a prime example. 



5 CONCLUSIONS 

We have analysed the high quality EPIC spectrum of NGC 45 17 
ULXl, demonstrating that this source exhibits spectral curvature 
at ~6 keV, similar to that seen in the highest quality observa- 
tions of other ULXs to date. This feature has been modelled with 
both the recent reflection and Co mp tonisation interpretatio ns (see 
ICaballero-Garcia & FabiaJlOlOl and lOladstone et al.ll2009l respec- 
tively) in an attempt to provide a direct comparison of the two. The 
reflection interpretation requires that the majority of the emission 
occurs within the innermost regions of the accretion disc. We find 
that the radial emissivity profile must be strongly centrally peaked 
in order to provide the smooth downturn, which can be achieved 
if the disc-coron a geometry is such that strong gravitational light 
bending occurs dMiniutti & F abian 2004). In order to model the 
curvature in NGC 45 17 ULXl, the disc must have, or appear to 
have, a highly super-solar iron abundance. This may be possible 
if the black hole is accreting from a debris disc formed primarily 
from the metal-rich inner layers of the star that originally collapsed 
to form it, or if the light bending is so strong that multiple reflec- 
tions of the intrinsic continuum occur. Meanwhile, we find that the 
Comptonisation interpretation requires that the electron tempera- 
ture of the Comptonising corona is very low, and that the corona is 
optically thick. It is the low temperature of the thermal component 
of the Comptonising electrons that provides the curvature in the 
spectrum, as this determines the high energy cut-off of the Comp- 
ton up-scattering process. These low temperatures and high optical 
depths may be due to mass loading of the corona i f disc winds are 
impor tant at high accretion rates, as proposed bv iKing & Pounds! 
( l2003h . 

Neither model may be excluded physically, so alternative ob- 
servational means are required to distinguish between them. Cor- 
rectly identifying the origins of the high energy curvature could be 
a very important step in enhancing our understanding of ULXs, as 
it is one of their characteristics that distinguishes them from XRBs. 
By extending both of the models to energies above 10 keV, it be- 
comes clear that there are large differences between them at these 
energies (Fig.|7]l. We expect this to be the case even if, as is likely, 
there is a contribution from a non-thermal distribution of electrons 
to limit the downward curvature in the Comptonisation model. The 
major cause of the discrepancy between the two models at high en- 
ergies is the presence of the Compton hump in reflection spectra, 
a broad emission feature that often peaks at ~30-50 keV. Obser- 
vations of ULXs above 10 keV should provide an important diag- 
nostic for the relative importance of these physical processes in the 
ULX population. 

Unfortunately, the angular resolution of the hard X-ray instru- 
ments aboard currently operational facilities that simultaneously 
provide spectral coverage over the ~0.5-10.0 and ~ 10-100 keV 
energy ranges is not sufficient to provide uncontaminated data for 
a large number of ULXs. We stress the importance of modelling 



data above and below 10 keV obtained simultaneously as spectral 
variability and accretion state transitions could otherwise lead to 
incorrect conclusions, especially if ULXs do primarily represent a 
new accretion state. Hard X-ray imaging spectrometers with good 
angular resolution, such as the Astro-H Hard X-ray Imaging Sys- 
tem, should allow a good number of reliable ULX spectra above 
10 keV to be obtained in the future. 
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